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26 Abstract 

27 Sexual antagonism occurs when there is a positive intersexual genetic correlation in 

28 trait expression but opposite fitness effects of the trait(s) in males and females. As 

29 such, it constrains the evolution of sexual dimorphism and may therefore have 

30 implications for adaptive evolution. There is currently considerable evidence for the 

3 1 existence of sexually antagonistic genetic variation in laboratory and natural 

32 populations, but how sexual antagonism interacts with other evolutionary phenomena 

33 is still poorly understood in many cases. Here we explore how self-fertilization and 

34 inbreeding affect the maintenance of polymorphism for sexually antagonistic loci. We 

35 expected a priori that selfing should reduce the region of polymorphism, since 

36 inbreeding reduces the frequency of heterozygotes and speeds fixation. Although this 

37 expectation was supported, our results show that there is an interactive effect between 

38 the degree of selfing and dominance such that those segregating sexually antagonistic 

39 loci that do exist are more likely to be partially dominant. In addition, inbreeding 

40 effects may influence population persistence and genomic location of sexually 

41 antagonistic loci in separate-sexed organisms. 
42 
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43 Introduction 

44 Sexual antagonism occurs when there is a positive intersexual genetic correlation in 

45 trait expression but opposite fitness effects of the trait(s) in males and females 

46 (Bonduriansky & Chenoweth, 2009; it is also known as intralocus sexual conflict 

47 when the intersexual genetic correlation is for the same trait in both sexes). As such, it 

48 constrains the evolution of sexual dimorphism and may therefore have implications 

49 for adaptive evolution. Sexual antagonism has often been considered a relatively 

50 transient phenomenon which will eventually be resolved by the evolution of sex- 

51 specific modifiers (Stewart et al, 2010), but recent research suggests that pleiotropic 

52 effects can constrain the evolution of sex-limitation (Mank & Ellegren, 2009), and 

53 that even the evolution of sexual dimorphism need not resolve the conflict (Cox & 

54 Calsbeek, 2009; Harano et al., 2010). Some authors have even suggested that sexual 

55 antagonism is inevitable in any species with separate sexes (Connallon & Clark, 

56 2014). The study of sexual antagonism has now matured to a point where its 

57 taxonomic ubiquity and potential importance in natural populations can no longer be 

58 questioned (Bonduriansky & Chenoweth, 2009; Cox & Calsbeek, 2009). The time has 

59 therefore come when we can no longer study sexual antagonism as an isolated 

60 phenomenon; it must now be integrated with broader evolutionary theory. To this end, 

61 we have constructed a model that investigates how self-fertilization and inbreeding 

62 affect the maintenance of polymorphism at sexually antagonistic loci. 

63 Our motivation was to begin formalizing recent speculations about the nature 

64 and relevance of sexually antagonistic variation in hermaphroditic species (Abbott, 

65 201 1; Bedhomme et al., 2009). All else being equal, sexually antagonistic alleles in 

66 hermaphrodites should exhibit the same dynamics as autosomal sexually antagonistic 

67 loci in separate-sexed organisms. This implies that selection must be strong and 
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68 approximately equal in magnitude across the sexes in order for polymorphism to be 

69 maintained (Kidwell et at, 1977). However hermaphrodites have a number of 

70 properties that set them apart from species with separate sexes, one of which is the 

7 1 fact that many hermaphrodites are partially or completely selfing (Goodwillie et al. , 

72 2005; Jarne & Auld, 2006). Selfing should reduce the region of polymorphism, since 

73 inbreeding reduces the frequency of heterozygotes and speeds fixation (Lynch & 

74 Walsh, 1998). One might therefore expect a priori that partially selfing 

75 hermaphrodites should exhibit lower levels of sexually antagonistic genetic variation 

76 compared to separate-sexed species or obligate outcrossers. However recent work by 

77 Fry (2009) demonstrates that dominance effects can have considerable influence on 

78 the region of parameter space permitting polymorphism. We therefore extended the 

79 framework developed by Kidwell et al. (1977) and Fry (2009) to investigate how 

80 proportion selfing and dominance interact to affect the maintenance of sexually 

81 antagonistic loci. We found that although inbreeding reduces the region of parameter 

82 space permitting polymorphism overall, it can offset some of the effects of dominance 

83 demonstrated by Fry (2009). 
84 

85 Model 

86 Our model is based on a classic framework for the investigation of sexually 

87 antagonistic alleles (Kidwell et al., 1977). Our population is made up of a large 

88 number of diploid hermaphroditic individuals. We focus on a single locus, at which 

89 there are two alleles, denoted A and a. This means that every individual is one of three 

90 genotypes: AA, Aa, or aa. Each genotype confers a different fitness; there are assumed 

91 to be no other differences between individuals. 
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92 We model generations as being discrete and non-overlapping. Within each 

93 generation, the life cycle goes as follows. We first census the genotypes in the 

94 population, and denote the frequency of genotype AA by p, and the frequency of 

95 genotype aa by q. Then the frequency of heterozygote Aa types is l—p — q. 

96 After censusing, random mating occurs. A proportion F of matings are self- 

97 fertilisation, while the remaining 1 — F matings are outbreeding events. For the self- 

98 fertilisation events there is no effect of genotype on offspring production. The 

99 genotypes of the offspring from self-fertilisation will depend on the parental 

100 genotype. Homozygous AA or aa individuals will produce their own genotypes for 

101 offspring, while heterozygous Aa individuals will have % of their offspring of 

102 genotype AA, Vz genotype Aa, and X A genotype aa (Table 1). Thus the frequency of 

103 each genotype in the next generation due to offspring from inbreeding events is as 

104 follows. AA:F(p + l A{\ -p- q)); Aa: Vi F(l -p - q); aa: F (q + Va{\ -p- qj). 

105 Because we are assuming random mating we can model each outbreeding 

106 event as being the result of the combination of sperm and eggs from randomly drawn 

107 individuals, with the probability of drawing a given genotype in a given sex role being 

108 proportional to the frequency of that genotype, and to its fitness in that sex role. 

109 Fitness differs across genotypes and across each sex role, as summarised in Table 2. 

1 10 The A allele is female -beneficial, male-deleterious, so that bearing an A allele makes a 

111 hermaphroditic individual better at the female role but worse at the male role. 

112 Conversely, the a allele is female-deleterious, male-beneficial, so that bearing an a 

113 allele makes an individual worse at the female role but better at the male role. These 

1 14 deleterious and beneficial effects are summarised by the parameters Sf and s m , which 

115 represent the selection coefficients, and /z/and h m , which represent the dominance 

1 16 coefficients (Table 2). 
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117 Since we assume that there are a large number of outbreeding events, the 

118 frequency of each genotype in the next generation due to offspring from inbreeding 

1 19 events is as follows. The frequency of AA individuals is 

Jp+y 2 {l P q)(l h fSf ))((l s m )p+y 2 {\ p q)(l h m s m )) 

120 (1 F)- — n — '-, 

f m 

121 the frequency of aa individuals is 
, P Q)(l h f s f ) + q(l s f ))(y 2 (l p q)(l h m s m ) + q) 



122 (1 F)- 



130 



f m 



123 where f and are respectively the mean fitness in the female and male roles, 

~ f = p+(l p q)(l h f s f )+q(l s f ), 

124 _ / \ 

m =P(l sj + (l p q)(l h f s f ) + q. 

125 The frequency of Aa individuals is the balancing expression so that these three 

126 frequencies add to (1 - F). 

Ill Putting together both self-fertilisation and outbreeding events, we can derive 

128 expressions for the change in frequency of genotypes AA and aa from one generation 

129 to the next, denoted Ap and Aq respectively, as 

P=F(p+/ 4 (l p q)) 
+ (l F) (p+^(1 P q)(l h f s f ))(p{l s m ) + y 2 (l p q)(l h m s m )) 

f m 

q=F(y 4 (l p q) + q) 
+ (l P # h f s f )+q(l s f ))(y 2 (l p q)(l h m s m ) + q) ^ 



f m 



131 (1) 

132 Using equations (1) we can establish whether the A and a alleles are protected from 

133 extinction when rare, and consequently whether polymorphism is protected or not 
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134 (Appendix), depending on the values of parameters F, Sf, s m , hf, and h m . When F = 0, 

135 we recover the classic results for this model (Fry, 2009; Kidwell et ah, 1977). 

136 Therefore in our analysis here we focus on the effect of F on the region admitting 

137 polymorphism. 
138 

139 Results 

140 We can use (1) to calculate expressions a and A corresponding to protection when 

141 rare of a and A, respectively. The female-deleterious, male-beneficial allele a is 

142 protected from extinction where rare if 

* ( F+ 2h '( l F )) S f 

143 s m >a = — \ — ; r/ \. 

2 F(l s f ) + 2(l F)(h fSf h m ) 

144 while the female-beneficial, male-deleterious allele A is protected from extinction 

145 when rare if 

, (F+2(l h,){l F))s f 

146 S„<A = 1 1 ; M 



(2^7 m +F(l 2h m ))(l S f )' 
147 We want to know the effect of F on these thresholds. Thus, we calculate 

a = 2(l h f h m )s f 

f (2 f(i Sf )+2(i F)(h fSf h m )y 

148 

A 2(h f + h m l)s f 



F (2h m+ F(l 2h m ))(l Sf ) 

149 If daldF < 0, and 8A /8F > 0, then increasing F increases the region of parameter 

150 space that leads to polymorphism. This is true if h/+ h m > 1, and false if hf + h m < 1. If 

151 hf+ h m = 1 (for example if both alleles are additive in their effects), changes in the 

152 proportion of inbreeding make no difference to whether or not there is a 

153 polymorphism. 
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154 Interestingly, this result relates to a previous finding, which showed that in the 

155 absence of inbreeding, the higher the value of h/+ h m , the smaller the region of 

156 parameter space admitting polymorphism (Fry, 2009). For any fixed value of F this 

157 remains the case in our model. However, this dominance effect is weakened by 

158 inbreeding, because inbreeding results in fewer heterozygotes, and consequently the 

159 effect of hf and h m is weakened (Figure 1). The results when F = 1 are identical to the 

160 case when hf = h m = 0.5. 

161 For weak selection (e.g. parameter values 0 < s/, s m < 0.1) there is very little 

162 scope for polymorphism when h/+ h m > 1 (Figure 2, see also Fry (2009)) regardless of 

163 the value of F. However, when h/+ h m < 1, the range of parameter values for which 

164 there can be a sexually antagonistic polymorphism due to weakly selected alleles is 

165 severely curtailed by self-fertilisation (Figure 2). 
166 

167 Discussion 

168 In hermaphrodites, the ability to self -fertilise will affect the maintenance or otherwise 

169 of sexually antagonistic polymorphisms. We expected a priori that selfing should 

170 reduce the region of polymorphism, since inbreeding reduces the frequency of 

171 heterozygotes and speeds fixation (Lynch & Walsh, 1998). Although this expectation 

172 was supported, our results show that there is an interactive effect between the degree 

173 of selfing and dominance. In species in which there is no self-fertilisation (e.g. 

174 separate-sexed organisms and obligate outcrossers) the more an allele that is 

175 deleterious to one sex is dominant in that sex, the smaller the region of parameter 

176 space that will admit polymorphism (Fry, 2009). However, this effect is weakened by 

177 self-fertilisation (Figure 1), so that in partially selfing hermaphrodites we would 

178 expect more dominant sexually antagonistic alleles remaining at polymorphism (and 
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179 fewer recessive alleles) than if there were no selfing. In particular, for weakly selected 

180 sexually antagonistic alleles which are on average partially recessive in their 

181 deleterious state, the range of parameter space allowing for a polymorphic equilibrium 

182 is strongly restricted in the case where there is inbreeding (Figure 2); if the alleles are 

183 on average partially dominant in their deleterious state, the region of parameter space 

184 allowing for polymorphism is increased by inbreeding, but remains small. Overall, 

185 therefore, the more self-fertilisation occurs in hermaphrodites, the fewer sexually 

186 antagonistic polymorphisms we would expect (assuming that most sexually 

187 antagonistic selection is weak). 

188 Our result is clear and simple, and shows how the ability to self-fertilise can 

189 have a strong effect on the genetics of a species. However, the effect of inbreeding is 

190 not strong enough to completely cancel out the effect of dominance: for a given 

191 frequency of inbreeding, it will still be the case that the more dominant the alleles are 

192 in their deleterious context, the smaller the region of parameter space in which they 

193 can exist at polymorphism. 

194 It is of course well-established that selfing can lead to inbreeding depression, 

195 but that in habitually selfing organisms the benefits of selfing should outweigh the 

196 costs of inbreeding depression (Goodwillie et ah, 2005). We have therefore not 

197 explicitly modelled inbreeding depression, and have assumed that selfed gametes do 

198 not experience selective effects of the sexually antagonistic alleles. The logic behind 

199 this assumption was that because selfing is usually considered a form of reproductive 

200 assurance (Goodwillie et al, 2005), this assurance will be ineffective if selfed 

201 gametes are subject to the same selection pressures as outcrossed gametes. This 

202 assumption is realistic if most of the cost of outcrossing is due to extrinsic factors, 

203 such as sexual conflicts with the mating partner (Anthes & Michiels, 2007; Koene, 
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204 2006; Koene et at, 2005) or energetic or predation costs of finding a mate (Jennions 

205 & Petrie, 1997). It becomes less realistic if the sexually antagonistic alleles cause 

206 intrinsic fitness differences (e.g. poor survival of gametes). Sperm (or pollen) 

207 limitation is unlikely to be a major limiting factor in fecundity when selfing (but see 

208 Hodgkin & Barnes, 1991), but it is not unlikely that mutations affecting egg 

209 quality/survival would have an effect even on the production of selfed offspring. If 

210 sperm are accompanied by toxic seminal fluid used in sperm competition when 

211 outcrossing, then this could also contribute to lower egg survival, even when selfing 

212 (Koene et ah, 2010; Scharer et ah, 2014). On the other hand, positive selection of 

213 alleles with a negative impact on egg production should be weak in species with high 

214 levels of selfing, because there is limited scope for fitness gain via sperm donation 

215 during outcrossed matings. We therefore believe that this assumption is reasonable. 

216 It is also worth noting that our model, although originally constructed with 

217 hermaphrodites in mind, is equally applicable to separate- sexed organisms with 

218 respect to inbreeding instead of selfing (Appendix). This generates some interesting 

219 predictions, especially for populations with high levels of inbreeding, such as island 

220 populations (e.g. Grant et ah, 2003), or populations with low dispersal levels due to 

221 habitat fragmentation (e.g. Andersen et ah, 2004). 

222 Firstly, species or populations with high levels of inbreeding should exhibit 

223 reduced levels of sexual antagonism compared to outbred populations, because the 

224 gain in polymorphism in dominant alleles is more than offset by the loss in 

225 polymorphism in recessive alleles, for a given level of inbreeding (Figure 1). A recent 

226 model suggests that sexual antagonism and demography can interact to cause 

227 extinction of populations located in patches that are beneficial to male fitness and 

228 detrimental to female fitness (Harts et ah, 2014). This is because populations collapse 
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229 if there are too few reproducing females. However populations which are declining in 

230 numbers should also become more inbred, as a result of the decreasing effective 

231 population size. Our results therefore suggest that the negative effect of sexually 

232 antagonistic alleles on population growth may be mitigated by inbreeding. This could 

233 reduce the chance of population collapse. 

234 Secondly, we should expect a higher proportion of sexually antagonistic 

235 alleles to be located on the X- or Z-chromosomes in inbred populations relative to 

236 outbred populations. Heteromorphic sex chromosomes have long been held to be 

237 likely hotspots for sexual antagonism. In a seminal paper, Rice (1984) argued that the 

238 X-chromosome should harbour increased levels of sexually antagonistic genetic 

239 variation because male-benefit loci that are recessive in females will be expressed in 

240 hemizygous males, but largely escape counter- selection in females at low to 

241 intermediate frequencies. Conversely, dominant female-benefit loci will also be more 

242 common on the X than on the autosomes, despite their deleterious effect in males, 

243 because the X spends more time in females than the autosomes (2/3 versus Vz) and 

244 therefore experiences stronger total female- specific selection. Our results suggest that 

245 the heteromorphic sex chromosomes should harbour higher levels of sexually 

246 antagonistic genetic variation in inbred populations, independent of the direction of 

247 the fitness effect (i.e. female-benefit/male-detriment versus male-benefit/female- 

248 detriment). This is because any X- or Z-linked locus that is not completely recessive 

249 in the homogametic sex will be partially dominant overall (i.e. h/+ h m > 1 will always 

250 hold true when homogametic > 0 because hheterogametic =1), and therefore subject to an 

251 increased range of polymorphism with increasing inbreeding level. 

252 In sum, we show that although inbreeding reduces the region of parameter 

253 space permitting polymorphism overall, it can offset some of the effects of dominance 
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254 demonstrated by Fry (2009). This means that although hermaphrodites with high 

255 levels of inbreeding are unlikely to harbour significant sexually antagonistic genetic 

256 variation, those segregating sexually antagonistic loci that do exist are more likely to 

257 be partially dominant. In addition, inbreeding effects may influence population 

258 persistence and genomic location of sexually antagonistic loci in separate- sexed 

259 organisms. 
260 

261 

262 Acknowledgements 

263 SJT was supported by funding from the European Research Council under the 7th 

264 Framework Programme of the European Commission (PBDR: Grant Agreement 

265 Number 268540), and JKA was supported by funding from the Swedish Research 

266 Council. 
267 

268 References 

269 

270 Abbott, J. K. (2011) Intra-locus sexual conflict and sexually antagonistic genetic 

271 variation in hermaphroditic animals. Proc. R. Soc. Lond. B Biol. Set, 278, 161-169. 
272 

273 Andersen, L. W., Fog, K., and Damgaard, C. (2004) Habitat fragmentation causes 

274 bottlenecks and inbreeding in the European tree frog (Hyla arbor ea). Proc. R. Soc. 

275 Lond. B Biol. Set, 271, 1293-1302. 
276 

277 Anthes, N. and Michiels, N. K. (2007) Precopulatory stabbing, hypodermic injections 

278 and unilateral copulations in a hermaphroditic sea slug. Biol. Lett., 3, 121-124. 
279 

280 Bedhomme, S. et al. (2009) How does breeding system variation modulate sexual 

281 antagonism? Biol. Lett., 5, 717-720. 
282 

283 Bonduriansky, R. and Chenoweth, S. F. (2009) Intralocus sexual conflict. Trends 

284 Ecol. Evol, 24, 280-288. 
285 

286 Connallon, T. and Clark, A. G. (2014) Evolutionary inevitability of sexual 

287 antagonism. Proc. R. Soc. Lond. B Biol. Set, 281, 20132123. 
288 



Downloaded from http://biorxiv.org/ on September 18, 2014 



289 Cox, R. M. and Calsbeek, R. (2009) Sexually antagonistic selection, sexual 

290 dimorphism, and the resolution of intralocus sexual conflict. Am. Nat., 173, 176-187. 
291 

292 Fry, J. D. (2009) The genomic location of sexually antagonistic variation: some 

293 cautionary comments. Evolution, 64, 1510-1516. 
294 

295 Goodwillie, C, Kalisz, S., and Eckert, C. G. (2005) The evolutionary enigma of 

296 mixed mating systems in plants: occurrence, theoretical explanations, and empirical 

297 evidence. Annu. Rev. Ecol. Evol. Syst., 36, 47-79. 
298 

299 Grant, P. R. et al. (2003) Inbreeding and interbreeding in Darwin's finches. Evolution, 

300 57,2911-2916. 
301 

302 Harano, T. et al. (2010) Intralocus sexual conflict unresolved by sex-limited trait 

303 expression. Current Biology, 20, 2036-2039. 
304 

305 Harts, A. M. F., Schwanz, L. E., and Kokko, H. (2014) Demography can favour 

306 female-advantageous alleles. Proc. R. Soc. Lond. B Biol. Set, 281, 20140005. 
307 

308 Hodgkin, J. and Barnes, T. M. (1991) More is not better: brood size and population 

309 growth in a self-fertilizing nematode. Proc. R. Soc. Lond. B Biol. Sci., 246, 19-24. 
310 

311 Jarne, P. and Auld, J. R. (2006) Animals mix it up too: the distribution of self- 

312 fertilization among hermaphroditic animals. Evolution, 60, 1816-1824. 
313 

314 Jennions, M. D. and Petrie, M. (1997) Variation in mate choice and mating 

315 preferences: a review of causes and consequences. Biol. Rev., 72, 283-327. 
316 

317 Kidwell, J. F. et al. (1977) Regions of stable equilibria for models of differential 

318 selection in the two sexes under random mating. Genetics, 85, 171-183. 
319 

320 Koene, J. M. (2006) Tales of two snails: sexual selection and sexual conflict in 

321 Lymnaea stagnalis and Helix aspersa. Int. Comp. Biol., 46, 419-429. 
322 

323 Koene, J. M., Pfortner, T., and Michiels, N. K. (2005) Piercing the partner's skin 

324 influences sperm uptake in the earthworm Lumbricus terrestris. Behav. Ecol. 

325 Sociobiol, 59, 243-249. 
326 

327 Koene, J. M. et al. (2010) Male accessory gland protein reduces egg laying in a 

328 simultaneous hermaphrodite. PLoS One, 5, elOl 17. 
329 

330 Lynch, M. and Walsh, B. (1998) Genetics and analysis of quantitative traits. Sinauer 

331 Associates, Inc., Sunderland, MA. 

332 Mank, J. E. and Ellegren, H. (2009) Are sex -biased genes more dispensable? Biol. 

333 Lett., 5, 409-412. 
334 

335 Rice, W. R. (1984) Sex chromosomes and the evolution of sexual dimorphism. 

336 Evolution, 38, 735-742. 
337 



Downloaded from http://biorxiv.org/ on September 18, 2014 



338 Scharer, L., Janicke, T., and Ramm, S. A. (2014) Sexual conflict in hermaphrodites. 

339 Rice, W. R. and Gavrilets, S. (eds). The genetics and biology of sexual conflict. Cold 

340 Spring Harbour Laboratory Press, Long Island, New York. Cold Spring Harbour 

341 Perspectives in Biology. 

342 Stewart, A. D., Pischedda, A., and Rice, W. R. (2010) Resolving intralocus sexual 

343 conflict: genetic mechanisms and time frame. Journal of Heredity, 101, S94-S99. 
344 

345 
346 



Downloaded from http://biorxiv.org/ on September 18, 2014 



347 Figure legends 

348 Figure 1 

349 The effect of inbreeding and dominance on the maintenance of sexually antagonistic 

350 polymorphisms in hermaphrodites. The area between two matching curves is where 

351 polymorphism is admitted. The green curves correspond to the case where h/=h m = 

352 0.1, so that the allele that is deleterious in each sex is partially recessive in that sex. 

353 The orange curves correspond to the case where h/= h m = 0.6, so that the allele that is 

354 deleterious in each sex is partially dominant in each sex. The dashed curves represent 

355 where F = 0, the situation where there is no inbreeding. The solid curves represent the 

356 case where F = 0.75, so that three quarters of matings are self-fertilisation. For the 

357 green curves, this results in a smaller area of polymorphism, while for the orange 

358 curves, it results in a larger area of polymorphism. The black dotted line is the 

359 asymptotic limit F = 1. It exactly corresponds to the case in which F = 0 and hf=h m = 

360 0.5. 
361 

362 Figure 2 

363 The effect of inbreeding on the maintenance of weakly sexually antagonistic 

364 polymorphisms. The region between any two matching-coloured lines admits a stable 

365 sexually antagonistic polymorphism. The pairs of matching lines correspond to the 

366 cases where F = 0, 0.25, 0.5, 0.75, and 1, respectively, as marked. Here h/= h m = 0.3; 

367 as inbreeding increases, the region admitting polymorphism decreases in size, to the 

368 limiting case where F = 1. For values of hf+ h m > 1, the region admitting 

369 polymorphism is contained within the region for F = 1, and consequently for these 

370 dominance parameters there is very little scope for polymorphism under weak 

371 selection. 
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372 Tables 
373 



Parental genotype 


Offspring genotype 


AA 


AA 


Aa 


l AAA, VzAa, l A aa 


aa 


aa 



314 Table 1: Offspring genotype distributions when inbreeding 

375 



Genotype 


Fitness in female role 


Fitness in male role 


AA 


1 


1 — Sm 


Aa 


!" V/ 


1 — h m S m 


aa 


l-Sf 


1 



376 Table 2: Fitness in different sex roles when outbreeding 

377 
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378 Appendix 

379 Stability of equilibria 

380 Using equations (1) we can define the function g[p, q] = (Ap, Aq), defined for all 

381 possible values of p and q (i.e. on the standard 2-simplex). We know that g[l, 0] = (0, 

382 0) (corresponding to fixation of the A allele), and g[0, 1] = (0, 0) (corresponding to 

383 fixation of the a allele). To determine whether either of these two equilibria are stable 

384 we consider the Jacobian matrix J of the function g, 

dp dp 

385 J: 



dp dq 
d q d q 
dp dq 

386 For each fixed point, we evaluate J and calculate its eigenvalues. If they are all 

387 negative for a given equilibrium point, that point is stable (thus if any of the 

388 eigenvalues are positive, the equilibrium point is unstable). If the equilibrium point at 

389 (1, 0) is unstable, then a is protected from extinction when rare (corresponding to the 

390 condition s m > a given in the main text). If the equilibrium point at (0, 1) is unstable, 

* 

391 then A is protected from extinction from rare (corresponding to the condition s m < A 

392 given in the main text). If both alleles are protected from extinction when they are 

393 rare, then we have a protected polymorphism. 
394 

395 Applicability of model to separate-sexed species 

396 Although the model was constructed to consider hermaphrodites, it can also apply to 

397 separate-sexed species. Because separate-sexed species cannot self-fertilise, the 

398 definition of F as the proportion of self-fertilising events cannot be maintained. 

399 Instead, F is taken to be a measure of the additional probability with which an 

400 individual will mate with a partner sharing the same genotype at the Ala locus of 
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401 interest (Appendix Table 1). Thus F can be seen as a measure of the level of 

402 inbreeding that is occurring in the population. 
403 



Focal 


Mating partner's 


Probability of that 


Probability of that 


genotype 


genotype 


partner genotype 


partner genotype with 






under random 


"self-fertilisation 






mating 


frequency" F 


AA 


AA 


P 


p + F(l-p) 




Aa 


1 -p-q 


(l-F)(l-p-q) 




aa 


q 


(1-F)q 


Aa 


AA 


p 


(X-F)P 




Aa 


l-p-q 


l-p-q + F(p + q) 




aa 


q 


(1-F)q 


aa 


AA 


p 


(l-F)p 




Aa 


l-p-q 


(l-F)(l-p-q) 




Aa 


q 


q + F(l-q) 



404 
405 



Appendix Table 1: App 



ication of F to separate-sexed species 



